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We extend the active walker model to address the formation of paths on gradients, 
which have been observed to have a zigzag form. Our extension includes a new rule which 
prohibits direct descent or ascent on steep inclines, simulating aversion to falling. Further 
augmentation of the model stops walkers from changing direction very rapidly as that 
would likely lead to a fall. The extended model predicts paths with qualitatively similar 
forms to the observed trails, but only if the terms suppressing sudden direction changes 
are included. The need to include terms into the model that stop rapid direction change 
when simulating mountain trails indicates that a similar rule should also be included in 
the standard active walker model. 
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1. Introduction 

The dynamics of pedestrians and their interactions with the environment have be- 
come a central theme in the study of social physics ®. An aspect of pedestrian dy- 
namics that has received considerable attention is the description of spontaneously 
formed trail systems. Unexpected patterns can be found in trail systems, highlight- 
ing the interplay between effective attractions and itinerancy when walking from a 
starting point to a destination R The creation of powerful models of human trail 
formation could help to improve planning and understanding of paths and related 
phenomena l 7 l 9 l 8 l 10 l 4 l 5 [ However, we are not aware of any physics based studies of 
the formation of mountain trail systems, where poor planning for human activity 
can lead to significant environmental damage 

In principle, there are a huge number of possible routes to be explored when 
choosing a path, and walkers could take any course between a starting point and 
a destination. A first approximation to the most probable path is a straight line 
between the initial and destination points (unless there are obstacles in the way). 
However, studies of trails using the active walker model have shown that the de- 
tailed patterns of paths form from a counterpoint between the desire to walk on 
well trodden paths, and the shortest route to be found by traveling directly be- 
tween the origin and destination EH Well-trodden paths are likely to be favored by 
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pedestrians because of reduced energy usage when compared, for example, to walk- 
ing through long grass. This preference may be largely psychological, as internet 
based experiments in a virtual environment have also shown that 'walkers' tend to 
favor well-used 'paths' S3 The preference to walk on regularly used paths leads to 
an effective interaction between past and present walkers, indicating that there is 
interesting physics involved in the formation of such trails. 

On inclines, there is a third influence. Walkers may ascend slopes diagonally if 
the gradient of the incline becomes too steep to ascend directly. Walking at an angle 
to the line of fastest ascent has the effect of decreasing the effective gradient of the 
incline, permitting travel up steeper slopes. On descent, it may not be possible to 
walk directly down a very steep slope without becoming unbalanced, again leading 
walkers to take a diagonal path. Our aim in this article is to determine if the 
walkers' desire to avoid steep gradients, in combination with the rules of the active 
walker model, can be used to simulate the zigzag paths that can be observed in 
mountainous regions. 

The energetics of walking on an incline have been discussed by Alexander in a 
simple model of bipedal locomotion ^ From energy considerations, walkers aim to 
change the angle of ascent from the vertical if a hill reaches a steep enough gradient. 
In Ref. 1 , human walkers are omniscient, and are able to assess the energetic outlay 
for an entire route. In this way, they can assess if more energy would be expended 
taking a shallower and longer route, or if a quick hike up a steep route would 
be more favorable. We consider that it is unlikely that the information about all 
possible routes is available for global decisions to influence route and planning, and 
expect that decisions are more likely to be local. Moreover, modeling has concluded 
that humans choose well trodden paths in a more local manner ®. We therefore 
consider that insight into mountain trail formation could be gained from active 
walker simulations of inclined planes. 

To demonstrate some examples of paths on inclines, we took photographs in the 
Lake District in Cumbria, England, which can be seen in Fig. [T| The upper panels 
Fig. [ljAl,A2) show spontaneously formed zig-zag paths on Wansfell. Two paths 
can be seen to the left and right hand sides of the picture, as highlighted in the 
figure on the right (A2). The left hand path has been augmented with rock since 
its formation, but still shows the characteristic zig-zag. To remove doubt on the 
origin of the zig-zags (for example, the stones might have been laid according to a 
plan) a second path can be seen to the right of the picture. The right hand path 
is spontaneously formed in the grass. Both have similar path angle and regularity 
of direction change. An example from further up the trail can be seen in panels 
Bl and B2. Here, there is additional wear to the side of the trail, indicating that 
the trail is still evolving. The paths are on inclines of around 1:8. The bottom path 
(CI and C2) is on an incline of around 1:2, and was observed near Low Sweden 
Bridge which is close to the village of Ambleside. Trails of this type are not unique 
to England, and such paths can be seen in other locations, such as Smith Rock in 
the USA, where the bends in the paths are large enough to appear on trail maps 
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Fig. 1. (Color online) Spontaneously formed zig-zag paths on Wansfell, near Ambleside, Lake 
District, Cumbria, UK. (Al, A2) The left hand path has been augmented by humans since its 
formation. The right hand path is spontaneously formed in the grass. Both paths have a similar 
angle relative to the steepest direction and regularity of direction change.). (Bl, B2) Another 
zigzag path from further up the trail. Recent wearing of the path can be seen. The top paths 
are on inclines of around 1:8. The bottom path (figures CI, C2) is on a much steeper incline of 
around 1:2. This path was found close to Low Sweden Bridge near Ambleside, and breaks off from 
the main path that can be seen curving around towards the back of the picture. White lines are 
included on the right hand versions of the figures to highlight the paths. 

m 

This article continues as follows. We review the rules of the active walker model 
in section [2] In section [3] we briefly discuss aspects of the biomechanics of walking 
on inclines. Extensions of the active walker model for mountain trail systems are 
introduced in section [4] Results of simulations are shown in section [5] Finally we 
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summarize in section [6] 

2. Active walker model for human trails 

In order to investigate trail formations on mountains we use a modified form of 
the active walker model, which was introduced in Refs. |7] and [5] In this section, 
we review the rules of the unmodified active walker model following the scheme 
introduced in Ref. [3 Since the active walker model forms the basis of our extension 
to mountain trails, it is our aim to ensure that all specific features of the active 
walker model are clear, before introducing our extensions to the model in section [4] 
The aim of the active walker model is to describe the formation of paths in soft 
ground. As pedestrians walk on soft surfaces such as grass, the ground becomes 
worn and a path emerges. In the active walker model, the wear on the soft surface 
is assumed to be represented by a function, G(r,<), which represents the ground 
condition at time t and position r. As is common in statistical physics, G(r,t) is 
assumed to evolve according to a first-order rate equation, 



The rate of change of G(r,t) depends on weathering (the first term on the right 
hand side of the equation) and wear by walkers (the second term). The rate of 
weathering is governed by the parameter 1 /T, where T sets the time scale of path 
decay back to the undisturbed ground condition Go- The saturation value is denoted 
G max . The parameter controlling the damage caused by a footprint is denoted /. 
Walkers are located at positions r«, and S(r) is the Dirac delta function. In the 
presence of more than one walker, \x is an index corresponding to the walker being 
considered. Note that G is a positive quantity in this study. 

The Dirac S- function found in equation [T] is not especially convenient to use for 
numerics. Therefore, we suggest the following small modification to equation [T] 



where f(r — rv) is an arbitrary function that specifies the shape of the damage 
caused by a footfall, /(r) is normalized to unity. We take / to be a square with 
sides I = 10cm long, i.e. to have similar area to the base of a shoe. With this choice 
of /, l 2 G max /I is the approximate number of footfalls that cause G to reach the 
saturation value G max . Since the feet of walkers have finite dimensions, we consider 
equation [2] to be more physical than equation [T] 

The positions of the walkers are updated according to how attracted they are 
to any local paths, and the direction of the destination. Each walker fj, is influenced 
by a potential that depends on their location relative to the paths. The attraction 
of local paths is determined from, 




(1) 
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The exponential function indicates how visible the path at r is to the pedestrian at 
location r^. The parameter, a controls the distance over which walkers are moti- 
vated to make detours to well formed trails. We note that the sign of the potential V 
is opposite to the usual convention, i.e. a positive V indicates an attraction. There 
are no negative V, since G is always positive. This is in agreement with the sign 
convention used in Ref. 9 

There are two competing influences on walkers in the active walker model. The 
first is the desire to reach the destination in as direct a manner as possible. In the 
conventional active walker model, each walker with index \x aims for a destination 
point, d M . In the absence of any paths, it is reasonable to expect walkers to travel by 
the most direct route, i.e. along the unit vector (d M — r M )/|d M — The attractive 
potential defined in equation [3] generates an effective force in the same direction as 
the gradient of the potential V r y(r^, t). In the active walker model, V r V(r(„ t) 
and (da — r„)/|d„ — rj are both used to define the direction e that walker [i moves 
along from position r M in the presence of a set of paths, 

, f) = (jg -r u (t))/K~ Ty(t)| + VrMr^t) 
e[T ^ > \(d,-r,(t))/\d,~r,(t)\ + V r ^(r,)\ [ > 

V r V(rn,t) is the direction of favorable ground. The symbol V r represents the 
gradient of V taken with respect to the position vector of the walker r M (the sub- 
script acts as a reminder that the gradient is not taken with respect to r). e is a 
dimcnsionless unit vector. 

The active walker model is completed with the equation of motion, 

^ = v 0li e(r^t) (5) 

where vq^ is the walking speed of walker /i, assumed to be constant for the entire 
journey of the walker. In the discretized form of the equations, vq^ is related to the 
length of a stride. 

The active walker model has been successful in describing some of the unex- 
pected features that have been observed in trail systems. The application of the 
active walker model to the winding paths found on steep inclines has not yet been 
considered. In the following section, we discuss the biomechanics of walking on in- 
clined planes with the aim of developing additional rules to explain the wiggles 
observed in trail systems on hills. 



3. The biomechanics of walking on inclines 

In this section, we aim to give a brief summary of the physiological factors that 
must be considered when modeling walkers on inclined planes. Some insight into 
the manner of walking that is adopted on an incline can be gained from attempts 
to develop bipedal robots that can climb sloping surfaces E^D. If exactly the same 
method is used to climb a slope as to walk on flat paths, the center of gravity is 
positioned further back than usual, which may lead to a fall^ (a good introduction 
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Fig. 2. Schematic of the region of discomfort for walking up or down the mountain. If the active 
walker model suggests angles in the shaded area \^y\ < 9, the angle will be mapped to either 8 or —8 
depending on which side of the axis the optimally desired angle lies. The mapping also depends on 
whether the walker is going up or down the incline. We estimate that a forbidden angle of around 
10° corresponds to a 1 : 4 gradient. Note that the x-axis is defined to be in the direction directly 
up the gradient, and the y axis along the contour of no height change. Angles are defined from the 
x-axis. 

to the biomechanics of walking on flat surfaces can be found in Ref . [T5|l . Humans 
(and robots developed to climb slopes) compensate for the effects of the slope by 
leaning forward, which can be achieved by flexing the ankle. This is satisfactory 
for very shallow slopes. However, if a slope is too steep, it becomes physiologically 
impossible to turn the ankle joint through the angle necessary to hold the center of 
gravity over the feet (for example, the person investigated in Ref. [TT|could not bend 
his/her ankle by more than 24'Q. Thus good foot contact cannot be maintained 
with the ground when walking directly up very steep inclines. 

Experiments have shown that the requirements on joints increase dramatically 
on moderate slopes l n l 12 l walking on a 10° incline requires hip flexibility of 60°, 
compared with 30° on the flat. Demands on ankle flexibility increase in a similar way. 
Joints in the leg such as the ankle, knee and hip are also subjected to significantly 
increased forces As determined in Ref. [TTJ the ankle becomes fully bent for more 
of the walking cycle as the gradient increases from 0° to 10°. The physiological 
constraints on the angles through which joints can bend indicate that hills eventually 
become too steep to walk up directly. To compensate for the limits of joint flexibility, 
the walker can choose to change the angle of ascent, to avoid walking directly uphill. 
This leads to a smaller effective gradient. The inability to walk directly uphill can 



a In biomechanics, the limit to bending of a joint is known as the maximum angular excursion 
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be thought of as an effective potential barrier against certain walking directions. 

There is also a risk of falling on descent. This is probably greater than the chance 
of falling on ascent, since parts of the walking cycle require the center of gravity 
to be moved forward ^ and on descent, the body is already unstable to falling 
forwards. Again, it is possible to compensate against falling by reducing the angle 
of descent. As before, we expect this compensation to lead to an effective potential 
barrier for angles that are likely to cause a fall (most walkers will have learned to 
avoid tumbling). When simulating mountain trails it is important to include this 
behavior. 

There is an additional feature of walking that should be relevant regardless of 
whether the walker is on an incline or level ground. Given the constraints of anatomy, 
sudden changes of direction between consecutive footsteps are likely to lead to 
discomfort, instability, falling or significant reduction in speed. Rapid changes of 
direction correspond to significant changes in the momentum of the walker and 
therefore large forces are required to maintain stability. These forces are applied 
at ground level, so can lead to a significant torque on the walker and toppling. 
We therefore expect walkers to have a preference to make consecutive footsteps in 
similar directions. We now describe the inclusion of incline effects and the inability 
to make sudden changes of direction into the active walker model. 

4. A model of mountain walkers 
4.1. New rules for mountain walking 

It is our aim extend the active walker model to construct a "mountain walker 
model" . As we have discussed in section [3] our model should take account of the 
inability of walkers on steep inclines to walk directly up or down a slope if the 
gradient becomes too great, and avoid sudden changes of direction that can cause 
instability. 

We suggest two ways of implementing an aversion to sudden changes in direction: 

Rule la The new direction of the walker is taken as the weighted average of the 
recent angle of motion <p and the angle of motion that would be favored on 
a flat surface. 

7 (t)=a0+(l-a)/3(t), (6) 

We call the parameter a the persistence of direction. [3 is the angle that is 
determined from the direction vector e = e x i + e y j, 

(3{t)= tan- 1 (e y (t)/e x (t)) (7) 

where i and j are unit vectors in the x (uphill) and y (horizontal) directions 
respectively. In the determination of /?, care is taken to ensure that the angle 
is in the quadrant consistent with the signs of e x and e y . We take care to 
avoid any problems with branch cuts in equations[6]and|7] Strictly, </> should 
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-2-10 1 2 



Horizontal [m] 

Fig. 3. Ground condition for (a) 9 = 15° (b) 6 = 25° and (c) 6 = 35°. In all cases, the persistence 
of direction a = 0. If consecutive steps are not taken in similar directions, no clear zig-zag paths 
form. 

relate to the direction of the walker over the time period just before 

<f>= [ (3(t')dt'/A (8) 

Jt-A 

where A is the time period within which the direction of previous footsteps 
is relevant to stability. 

We also suggest an alternative to equation [6j 
Rule lb Equation [2] is supplemented with the average direction over recent time, 
= K ~ *y(t))/|d M - r„(t)| + Vr.VM + 1 it a t')df 
IK " r„(t))/|d^ - r„(t)| + VrMrJ + f //_ A e(r M ,t')di'| ' 

(9) 

where again A is the time within which sudden change of direction would 
lead to falling and a is a parameter that controls the influence of the pre- 
vious footsteps. 

We also modify the active walker model to incorporate the potential barrier in 
angle space that represents the inability to achieve a safe journey if the walker is 
traveling too close to the directly uphill or downhill routes (i.e. a journey where the 
walker is unlikely to fall over) . This rule is summarized in figure y2\ 
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-2-10 1 2 

Horizontal [m] 

Fig. 4. Ground condition for simulations with minimum safe angle 8 = 5° and persistence of 
direction (a) a = 0.4 (b) a = 0.5 and (c) a = 0.6. As the tendency to make consecutive footfalls 
in the same direction increases, the size of the features increases. 

Rule Ila (To be used with rule la) If the angle 7, chosen from equation [6] falls 
within a forbidden zone (I7I < 9), the angle along which the walker moves is 
mapped to the nearest minimum angle of safety (either 9 or —6). We assume 
that the larger the incline of the slope, the larger the region of angles that 
are not permitted. This is equivalent to having an infinite effective potential 
barrier in certain directions. In the event that 7 = our walker has a 
preference to move left. If there is a mapping, we reconstruct e from 7 as 
e = cos(7)i + sin(7)j, noting that the sense of i is up the plane. The new 
e is then used to update the position of the walker via eqn. [5j 

Rule lib (To be used with rule lb) If the angle, (3 = tan -1 (e y j e x ) , lies within a 
forbidden zone (\f3\ < 9), then the vector e is mapped to the closest angle 
outside the forbidden region. If (3 = the walker moves left. Since j3 — 
on the first iteration, our paths show a bias in that direction. 

4.2. Discretization scheme 

We use the following simple first-order discretization procedure to solve the active 
walker model and extensions. The time coordinate is discretized into small intervals 
At, and the spatial coordinates into intervals Ax and Ay. Therefore, the vector r 
can be written in the discrete form r = (r/Ax, vAy) where 77 and v are integers. 
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Horizontal [m] 



Fig. 5. Ground condition for simulations with 9 = 15° and (a) a = 0.4 (b) a = 0.45 and (c) 
a = 0.5. 



Approximated using this scheme, equation [2] becomes, 
G(r,t n ) := G(r,t n _!) + ^ [G - G(r,t n ^)} 



■Y, I At 



1 - 



G(r,t n _i) 

G max 



/(r-r M (t n _i)), (10) 



The integral in equation [3] becomes a sum, so the potential can be computed as, 
V„(r,t n ) := ^AxAye-l Aa; '' i + A ^-^( t "- 1 )l/ ff G(Ax7 7 i + Ayvj,t n ), (11) 



77,1/ 



From the potential, the direction of the walker e can be computed. For convenience, 
we will write the numerator of eqn. [4] as the un-normalized vector e, and then nor- 
malize it once it has been computed to calculate e. The un-normalized x-component 
of eqn. [4] becomes, 



(d x -r x ) [V^x^ + Ax, Up) - V^x^ - Ax,y p )] 



\d-r\ 2Ax 
and the unnormalized y-component of eqn. [4] is approximately, 

( d y- r v) , [Vn(x li ,y fl + Ay) - V^x^y^ - Ay)] 



Syfj,(t n ) 



\d-r\ 



2Ay 



(12) 



(13) 
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Horizontal [m] 

Fig. 6. Ground condition for simulations with 9 = 25° and (a) a = 0.4 (b) a = 0.45 and (c) 
a = 0.5. 

Here and j/ M are defined via r M = x^i + y^j. The components e xfi and e yfl are 
combined into a single vector e = e xll i + e yfJi j which is then normalized to determine 
the direction of the next step, 

e(t n ):=e(t n )/\e(t n )\ (14) 
The average angle (equation [8]) is approximately, 

A/At 

<f> = P(tn-i)M/A (15) 

with a similar scheme for the average of previous directions in equation [9] Equation 
[9] is otherwise discretized in the same way as equation |12| 

The walker position is updated according to the discretized version of eqn. [5j 

r>(i n ) := r„(t n -i) + Atvoe^tn). (16) 

The vectors r „ were permitted to take any real value to avoid truncation errors 
that could lead to limit cycles. 

It is important to note that, in general, the validity of first order integration 
schemes is not guaranteed. We checked our code and approximate scheme against 
results in Refs. [7] and GO finding no significant differences. We have run simulations 
using 3 different resolutions of Ax, Ay and At (see the results section of this article). 
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The scheme was found to converge quickly on reduction of the size of the discrete 
steps. 

As humans tend to walk with different step sizes, we varied the speed of the 
individual walkers. This variation leads to continuous paths. The width of the area 
studied was chosen as 10m and the length of the observed area (in the direction of 
the incline) as 25m. This is probably realistic for this type of trail formation, since 
we suspect that walkers aim for local goals, rather than the final goal of the peak 
(which is not necessarily visible from all locations) . 



4.3. Algorithm one 

To clarify our computational scheme further, we detail the algorithms that we use 
for computation. Two computational schemes have been used. Algorithm one is 
based on Rules la and Ha, with algorithm two based on Rules lb and lib. The steps 
in the algorithm are: 

(1) To initialize, set the ground condition to the default value G(r, t ) = 0. 

(2) Choose a walker speed following equation = (0-5 + r)ms -1 , where r is a 
random number between and 1. 

(3) Select the direction for the walker (i.e. up or down-hill) randomly with equal 
probability. 

(4) Initialize the walker position at the bottom or top of the incline depending on 
whether an up or down walker has been selected. 



(5) Compute the attractive potential V from equation 11 

(6) Determine the default walker direction from equation |14| 

(7) Apply the persistence of direction formula, equation^ 

(8) Apply the forbidden angle rule. If the angle, 7, lies within a forbidden zone 
(I7I < 9), then the vector e is mapped to the closest angle outside the forbidden 
region. In the event that 7 = our walker has a preference to move left. 



10 



(9) Update the walker position using equation 16 

(10) Update the ground condition using equation 

(11) Repeat from step [5] until the walker has reached the destination. 

(12) Repeat from steppluntil the action of sufficient walkers on the ground condition 
has been investigated. 



4.4. Algorithm two 

(1) Initialize the ground condition. 

(2) Select the speed of the walker. 

(3) Select the direction of the walker. 

(4) Initialize the walker position. 

(5) Compute V from equation 11 

(6) Apply the alternative formula for persistence of direction, equation [9] 

(7) Apply the forbidden angle rule. If the angle, /?, lies within a forbidden zone 
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(\(3\ < 6), then the vector e is mapped to the closest angle outside the forbidden 
region. In the event that 7 = our walker has a preference to move left. 



(8) Update the walker position r M using equation 16 



(9) Update G using equation 10 



(10) Repeat from step p\ until the destination has been reached. 

(11) Repeat from step until the desired number of walkers have traversed the 
slope. 

The simulation was implemented in C++. The random number generator, ran2 
from 'Numerical Recipes in C was used to ensure a long period^. We tested our 
code against the triangular and square examples in Ref. [5] 



5. Results 

5.1. Algorithm one 

In this section, results from the mountain walker extension to the active walker 
model are shown. We begin by considering only walkers traveling down the slope. 
The following parameters were used for the simulations: the maximum ground po- 
tential was set to G max = 200m~ 1 (the units of G are set by Eq. [3] and Eq. [5]) while 
the minimum ground potential was set to Go = 0. Larger values of G max lead to 
larger attraction to the path. A new walker descended the inclined area as soon 
as the existing walker reached the destination (tests showed that there was little 
difference between this scheme and starting a walker every 100s). The visibility, a 
of the paths was set to 10m and the intensity, / was set to / 2 G max /./V, where N is 
the number of footprints needed to wear the ground condition to 1/e of its max- 
imum value. N was set to 50 footfalls and the weathering parameter was initially 
set to T = 1000s. The values of TV and T are smaller than in a real trail system, 
where N would be of the order of several hundred footfalls and T would be of the 
order of a few days. This still leads to realistic simulations since Helbing et al. have 
found that a combination of several parameters of the active walker model could 
be represented by the single parameter k = IT j a = G max T/N<T 7 . Individual walk- 
ers were assigned a random speed between 0.5m/s and 1.5m/s. In the simulation, 
walkers were given the starting position ri n ;ti a i = (0m, 5m) and a destination of 
'"final — (25m, 5m) with the x-direction being the distance down the incline and 
y-direction the distance across the incline. Initially, 25000 walkers traversed the 
incline in each simulation. 

In our initial simulations, we set the parameter a = 0, consistent with the 
standard setup of the active walker model. This led to a curious result, as shown 
in Fig. [3] We found that walkers changed their walking direction frequently, with 
some walkers changing direction on every time step. This waddling gait is not very 
comfortable on an incline, and it is likely to make the walker very unstable. Also, the 
speed of advance down a hill is likely to be quite slow. We thoroughly investigated 
the parameter space to look for the zigzag paths, but did not find them. As this 
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waddling gait is not a practical way to descend a mountain, it is reasonable to 
conclude that the active walker model is missing some crucial physics: a walker 
is much more likely to ambulate by moving consecutive legs in similar directions, 
otherwise the resulting motion will be very unstable. For this reason, we introduce 
a factor, a, to represent persistence of direction (see Eq. ^J. a controls a weighted 
average between the angles on consecutive time steps. 

By changing the weighting factor, which alters the influence of the recent direc- 
tion of motion on the walker, and by changing the size of the forbidden angle 9, we 
were able to estimate how large a needs to be to produce convincing zigzag paths. 
Again, we assumed that walkers always traveled down the incline (in the direction 
of increasing x) and otherwise used the same parameters. We ran simulations with 
9 ranging from 9 = 5° to 9 = 45° and a ranging from a — 0.1 to a — 0.9. 

We show the results of single runs for a minimum safe angle of 9 = 5° in figure [4] 
The zig-zag paths are immediately apparent. Multiple turns in the path can be seen 
when the parameter, a > 0.5 (Fig. |3Jb)). This indicates that when walkers choose 
where to place the next footfall, around half of the choice comes from a desire 
to maintain the previous walking direction, and half from the wish to reach the 
destination quickly via a path. The effect of reducing the influence of the previous 
angle can be seen in Fig. |4ja) where a — 0.4. As a result of the reduced influence 
of the previous direction, the walker's motion follows only a slight wiggle. In panel 
(c) of Fig. [4] where a — 0.6 and the influence of the previous angle is increased, the 
overall size of features increases. This is a common feature of all simulations where 
the influence of previous step directions is increased through a. 

As the minimum safe angle is increased to 15°, we note that the minimum values 
of a which are needed to generate zig-zag patterns decreases slightly. In Fig. [5ja) 
where a — 0.4, we observe that a single zigzag emerges in the path. As walkers get 
closer to their destination, the rate at which they change their walking direction 
increases. Fig. |5jb) shows a similar path for a — 0.45. In Fig. [5^c) we observe 
that when a — 0.5 walkers take larger detours away from the direct path between 
entry and exit, with several changes of direction along the trail. Similar results were 
found when 9 = 25° (Fig. [6|. While zigzag forms develop during the simulations, 
the rather diffuse paths are not satisfactory representations of mountain trails. 

We complete our simulations using algorithm one, by examining walkers that 
travel consecutively up and down the incline. Results of the simulations can be seen 
in figure [7] A slightly different weathering time of T = 1500s was used. There is 
no requirement that the forbidden angles are the same for hikers moving up and 
down the incline. In general, we expect that the different mechanisms for walking 
up and down inclines lead to different forbidden angles. In this set of simulations, 
we choose a forbidden angle of 9\ — 25° for walkers moving up the diagram (down 
the gradient). Those moving down the diagram (up the gradient) have a minimum 
safe angle of #2 = 10°. We run simulations for a range of a. Some zig-zag patterns 
can be seen for low alpha, but patterns of decent size do not appear until a > 0.45 
(panel d), Again, the sizes of the bends in the path increase in size with a. The 
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Fig. 7. The results of simulations of 2500 walkers traveling consecutively up and down the incline. 
Algorithm one is used (a) a = 0.30 (b) a = 0.35 (c) a = 0.40 (d) a = 0.45 (e) a = 0.50 (f) a = 0.55 
(g) a = 0.60 (h) a = 0.65. T = 1500s, N = 50, a = 10m and v = (0.5 + r)ms- 1 where r G [0, 1) is a 
random variate. We also considered simulations with 25000 walkers, with no noticeable difference 
to the pattern of paths formed. At = 1.0s 
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Fig. 8. As Fig. [7] with double the space and time resolution used in the computation Ax = Ay 
5cm and At = 0.5s. The shapes of the paths are robust against the finite size scaling. 
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combined effect of walkers moving in both directions are trails with zig-zag patterns 
that have a similar angle to the largest of the two forbidden angles. The walkers 
with a smaller forbidden angle round off the sharp turns in the paths that were 
found on the zig-zags formed when walkers are only permitted to move in a single 
direction. This rounding is a direct consequence of the attraction term in the active 
walker model and may explain the curved nature of spontaneously formed mountain 
trails. Moreover, the inclusion of walkers with different minimum safe angles leads to 
well defined trails (rather than the diffuse trails found previously) . This is probably 
because the walkers with smaller forbidden angles have more freedom to change 
direction, allowing the active walker rules to function effectively. 

To test the effects of discretization, we compute results using smaller At = 
0.5s, Ax = 5cm and Ay = 5cm in figure [8] In figure [9] we use even smaller steps 
At = 0.25s, Ax — 2.5cm and Ay — 2.5cm. The absence of new features shows our 
calculations to be stable against finite size scaling. 



5.2. Algorithm two 



We also run simulations using algorithm two. The results are shown in Figs. 10 



11 and 12 with the resolution of the discretization increased between figures. 2500 



walkers traversed the incline in both directions. The parameters of the runs were 
T = 1500s, N = 50, a = 10m and v ail = (0.5 + r)ms" 1 where r e [0, 1) is a 
random variate. a was varied. Again, we choose 9\ — 25° and 9-2 = 10°. The results 
from both algorithms one and two are in good qualitative agreement. Note that the 
results for a — 2.5 required slightly higher resolution for convergence. In comparison 
to the patterns from algorithm one, the zigzags are quite regular. Small wiggles of 
around 0.5metres across form when a = 1.5. As before, the features become larger 
as a is increased. 



6. Summary 

We have developed an extension to the active walker model to handle the formation 
of trails on inclines. Our simulations are in qualitative agreement with empirical 
observations of mountain path formations. Such trails are characterized by a zig- 
zag pattern. Our extension took account of the inability of walkers to walk directly 
up or down very steep gradients. We have shown that some supplementary rules 
need to be included in the active walker model to achieve features consistent with 
mountain trails. Those additional rules are that (a) the consecutive steps of walkers 
tend to be in the same direction and that (b) there is a maximum permitted angle of 
motion down an incline to avoid falling and (c) there is a maximum angle of ascent 
for physiological reasons such as limited ankle flexibility. When rules to encourage 
consecutive steps are not present, we find that walkers travel in an unusual manner, 
only taking a single step before changing direction. We also found that the presence 
of walkers moving both uphill and downhill (with different forbidden angles) is 
important for forming well defined zig-zag paths. Walkers traveling downhill with 
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a larger forbidden angle are constrained to form zigzags but create rather diffuse 
paths unless those paths are complemented by walkers traveling uphill. This happens 
because the walkers with more angular freedom become attracted to and reinforce 
the paths. 

The preference of walkers to take consecutive steps in the same direction could 
also be relevant for the understanding of trail formation in flat areas, and it is 
our opinion that a persistence of direction parameter should also be included when 
simulating trail systems on level ground. Thus, our main conclusions are that (1) 
correlation between consecutive step directions should be included in the active 
walker model and that (2) an extension of the active walker model to include for- 
bidden angles, correlation between consecutive steps and the combination of walkers 
moving both up- and down-hill is suitable to understand the formation of mountain 
paths. 
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Fig. 9. As Fig. [7] with four times the space and time resolution used in the computation Ax = 
Ay = 2.5cm and At = 0.25s. Again, the shapes of the paths are robust against the finite size 
scaling. 
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Fig. 10. The alternative direction averaging scheme for the mountain walker, (a) a = 0.5, (b) 
a = 1.0, (c) a = 1.5, (d) a = 2.0, (e) a = 2.5, (f) a = 3.0. At = Is, Ax = Ay = 10cm. Results 
from both schemes are in qualitative agreement. The size of features increases monotonically with 
a. 
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